Tunneling with dissipation and decoherence for a large spin. 
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We present rigorous solution of problems of tunneling with dissipation and decoherence for a spin 
of an atom or a molecule in an isotropic solid matrix. Our approach is based upon switching to a 
rotating coordinate system coupled to the local crystal field. We show that the spin of a molecule 
can be used in a qubit only if the molecule is strongly coupled with its atomic environment. This 
condition is a consequence of the conservation of the total angular momentum (spin -I- matrix), that 
has been largely ignored in previous studies of spin tunneling. 

PACS numbers: 75.45. +j, 75.80.-(-q 
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The problem of tunneling of a large 
spin lll|,||,||J|J|,§,0,||j|] has received considerable attention 
lately in connection with spin-10 magnetic molecules, 
Mn-12 and Fe-8, 0,|ll|,0,||jl|Jl|,|l|,|l3 . High-spin 
molecules have been proposed as qubits for quantum 
computers [p8| . It is, therefore, important to understand 
the effect of the atomic environment on spin tunneling 
and decoherence. In this Letter we present rigorous so- 
lution of both problems for a spin of a molecule in an 
isotropic solid. 

The Caldeira-Leggett approach jl^] to the problem of 
spin tunneling with dissipation due to phonons was out- 
lined in the second of Refs. It has been applied in 
Ref. , though a general expression for the effective ac- 
tion has not been obtained. The correct formulation of 
the problem should account for the conservation of the 
total angular momentum: spin S plus the angular mo- 
mentum L of the atomic lattice |2^. In the absence of 
the external field, tunneling of S between fl" and -IJ. should 
be accompanied by a simultaneous co-flipping of L and 
is possible only if J = S + L = 0. This can be satisfied by 
an integer but not by a half-integer S, which is another 
way to look at the Kramers theorem. The necessity to 
have L — S in the tunneling state results in the mechani- 
cal rotational energy Er — {TiSY /2X where I ^ pl^ is the 
moment of inertia of the atomic lattice; p and I being the 
mass density and the linear dimension of the solid matrix 
containing S. This energy should be compared with the 
tunneling splitting A. The condition A > Er is needed 
for the tunneling state with L = 5, J = to be the 
ground state of the system, while in the opposite case of 
A < Er the ground state should he L = 0, J = S, with 
S frozen along the anisotropy axis. This translates into 
a minimal size / (typically of order 1-10 nanometers) of a 
free particle whose spin can tunnel between equilibrium 
orientations. 

The problem of the decoherence is even more subtle. 
On one hand, the matrix elements of the conventional 



spin-lattice interaction due to the electrostatic crystal 
field vanish between (|_fr> +| and (| fr> -| 4J.>) 
tunneling spin states [p2| . On the other hand, the real- 
time coherent oscillations of S in a solid matrix must 
be accompanied, through the conservation of the angu- 
lar momentum, by the oscillating shear deformation of 
the solid. This should be the source of decoherence. 
The relevant size of the solid involved in that process is 
Ic = Ct/uJc, where Ct is the velocity of the transverse sound 
and ujc = A/h is the oscillation frequency. We shall see 
that r = Er/h^hS^/pl^ emerges as a frequency scale 
that determines the decoherence rate. 

The imaginary-time action of the system consists of 
the action of the spin in a crystal field of the atomic 
lattice and the action of the lattice. The spin action can 
be written as 



/,[n] ^Iwz[n]+ / dTE,{n) , 
Jo 



(1) 



where Iwz is the Wess-Zumino action Es is the en- 
ergy of the magnetic anisotropy due to the crystal field, 
n(r) — S{t)/S, and /3 — 1/T. The anisotropy is deter- 
mined by the local atomic environment of the atom or 
molecule with spin S. The global symmetry of the lat- 
tice is less important since the wavelengths involved in 
the problem (see below) are always large compared to 
the atomic scale. For that reason, and in order to sim- 
plify calculations, it is convenient to choose an isotropic 
(e.g., amorphous) solid that is characterized by two elas- 
tic moduli only, fi and A, which determine velocities of 
the transverse and longitudinal sound, ct — (fi/p)^^^, 
ci = [(A -I- 2/i)/p]^/^. Within linear elastic theory the 
corresponding phonon action is 



I,\u] = 



hp 



dr / d-^v\ -pii^ + ^{u.^Y + 



A 



(2) 



where u(r) is the phonon displacement field and Ui 



is the strain tensor. 
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We shall assume throughout this Letter that the in- 
teraction of the spin with its atomic environment is con- 
tained entirely in the structure of the crystal field re- 
sponsible for the magnetic anisotropy. For example, in 
Mn-12 the spin Hamiltonian is dominated by the uniax- 
ial anisotropy, Hs = —D{e-S)'^ where e is the unit vector 
along the anisotropy axes. The conventional way to intro- 
duce the spin-phonon interaction |jll|,|l^ is to consider a 
small perturbation of e by phonons, 6e = [S(j)xe], where 
S(f) — iVxu is the local rotation. This results in the 
spin-phonon Hamiltonian 

Hsp^ D{Sz,Sx}uJzx + D{Sz,Sy}uJzy , (3) 

where {Si,Sj} is the anti-commutator and LUij = 
^{diUj — djUi). The effective spin action, 

lef f [n] = [n] + lenv [n] , (4) 

should be obtained by computing the environmental ac- 
tion lenv, 



exp 



where 



Ienvm\ 1 ( /mt n, Ul 

= — (p Puexp 



ZtJ 



Zi= <p X>u(r, t) exp - 



n 



(5) 
(6) 



The conventional approach [^,^ is to choose lint = 
drl-Lspin) with Hsp of Eq.(3) or similar. Then the 
integration in Eq.(5), though cumbersome, is Gaussian 
and, in principle, can be performed exactly. As has been 
discussed above, the problem with such an approach is 
that the tunneling of S from the it state at r = — oo to the 
4 state at T = +00 must formally involve the co-flipping 
of L, that is, the change of the mechanical rotation of 
the solid from clockwise to counterclockwise. In a fixed 
coordinate system the corresponding instanton involves 
large dispacements u which are difficult to work with. 

The difficulty mentioned above can be avoided by 
switching to a coordinate system that is centered at the 
spin and is firmly coupled to the local anisotropy axes. 
It rotates in the presence of the time-dependent shear 
deformation. In such a coordinate system the magnetic 
anisotropy Eg [n] remains unaffected by phonons and the 
tunneling of the spin is accompanied by small lattice dis- 
placements only. The new term appears in the energy, 
though. 



E' = -ns-fi 



(7) 



where fl — S<j> = iVxii. This can be considered as the 
consequence of the fact that rotation is equivalent to the 
magnetic field. In the rotating coordinate system in 
Eq.(5) becomes 



— —ihS 



hp 



drh 



ivxu 



(8) 



where we have integrated by parts to move the time 
derivative to n. Consequently, Eq.(5) is again a simple 
Gaussian integral on phonon variables. 

Equations (7) and (8) can be mistakenly taken for a 
parameter-free spin-phonon interaction pl . It is impor- 
tant to understand, therefore, that they are solely the 
consequence of switching to a rotating coordinate sys- 
tem, where the interaction between S and the lattice is 
contained, to all orders on u, in the anisotropy energy 
i?s[n]. To illustrate the equivalence of this approach to 
the conventional method of studies of the spin-phonon 
interaction, let us compute, e.g., the width, Fi, of the 
T = spin-precession resonance for Hg = —DS^, that is 
the rate of the one-phonon decay \m ^ S — 1 >— > |m = 
S >, provided by Eq.(3) and by Eq.(7). After stan- 
dard quantization of u (see below for details) Eq. (3) gives 
Fi = 5(25'- l)^D^ujf/12Trhpcl, while Eq.(7) 
gives r[ — Shuj I / 1271 pc^, where huji is the distance be- 
tween the m — S and ni — S — 1 levels. Observing that 
this distance equals (25 — 1)D, one immediately obtains 
Fi = T'^. The equivalence of the two methods in gen- 
eral can be traced to the equation of motion for the spin 
operator, h{dS/dT) — [H,S]. 

We shall now proceed to the computation of the 
Caldeira-Leggett action. In terms of phonon modes. 



u(r,T) = ^^e*'' ''ekAWk,A(r,T) , 

^ kA 



(9) 



/; and lint become 



Il[u\= I dr ^M(M_kAMkA + ^kAM-k.A^kA) 
kA 

/i„t[u, n] = -i dr^h-CkAWkA , (10) 
•^0 n 



where 



CkA = — 7=?i5(ik)xekA 



(11) 



ek,A are phonon polarization vectors, e_k,A'ekA' = Sw , 
M is the unit cell mass, and N is the number of unit cells 
in the lattice. The result of the integration over phonon 
variables in Eq.(5) can be presented as a sum of the local 
and nonlocal terms. 



Ienv[n] = linvln] + lenvi^] , 

where the local part is given by 

n = — -— / drn , 



l2pVc Jo 
and the nonlocal part is 



(12) 



(13) 



Tnl 



32nficf Jq 



dr 



dr'nu^r-r T"'' ■ (w) 
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Here Vc — M/p is the unit cell volume, lj^ = ct/Vc^'^ is 
the Debye frequency for the transverse phonon branch, 
and 



,3 - 

X e 



(15) 



The typical scale of the time derivative in Eqs. (13) 
and (14) is set by the instanton frequency, uji, which also 
determines the temperature, Tc'^huji, of the crossover 
from superparamagnetism to quantum tunneling of S. 
The scale of Tc is set by iS^/S* ||^. It is easy to see that 
l^mv/^invl ^ {Tc/TdY- This ratio is small as long as the 
energy of the magnetic anisotropy is small compared to 
the elastic energy, which is normally the case. The rel- 
ative effect of the atomic environment on spin tunneling 
is given by the ratio lenv/Is, which is of the order of 
Es/fj,Vc. For a molecule of spin S, rigidly imbedded in a 
solid matrix, this ratio is typically 10~^ — 10"'^. On the 
contrary, a loose coupling of the magnetic molecule with 
its atomic environment is equivalent to — > 0, which must 
result in a drastic increase of lef / and the corresponding 
exponential decrease of the tunneling rate. 

The advantage of our method becomes apparent when 
one wants to compute the rate of the decoherence for 
quantum oscillations of S between ff and -IJ.. Here we 
assume that the tunneling splitting A has been already 
renormalized by the effects studied above. If A is small, 
the problem can be trancated to the spin 1/2 problem 
with only two states: the ground state |0 >= "^(Itt > 
+ I ^ >) and the excited state, |1 >= ^(| f|- > - | -IJ. >). 
At T = the decoherence occurs due to the transition 
|1 > 10 > via spontaneous emission of a phonon of 
energy A — huc- Such transitions seem to be absent 
for the conventional magneto-elastic coupling that one 
obtains by considering the adiabatic rotation of the elec- 
trostatic crystal field (see, e.g., Eq.(3)). To compute the 
decoherence rate, one has to take into account the dy- 
namical nature of the phonon [ p^ . This can be done by 
switching to a rotating coordinate system coupled to the 
local crystal field. 

The general formula for the transition rate is 



T = —J2<^\V\J><J\V\^>SiE, 



(16) 



Substituting here V — — TlS-fl, we get 

r = 2;i< o|s|i > [j(A)] < i|s|o > , (i7) 

where J(A) is the spectral function of the environmental 
coupling, 

J(A) = < kA|f2|0 >< 0|f2|kA > (5(A - hu^^^) . 



kA 



(18) 



r = 2hs^j{A) . (19) 

To compute the spectral function, let us consider 

J{A) =TrY,^lx'»^i^xH^~^^kx) , (20) 



kA 



where 



Here 



Qkx =< o\n\kx > 



2 Zp 



(21) 



(22) 



and n(r, t) is the momentum of phonons in terms of the 
operators of creation and annihilation, 




[akAekAe^"-- - h.c] . (23) 



This gives 

J^kA = 



TltJkA 1 



ikxBkx ■ (24) 



Consequently, Eq.(20) becomes 

J{A) = ^Y. ^^kA® CkA<^(A - hu;^x) ■ (25) 



kA 

The spectral function is J = ^Tr[J'{A)], which gives 

'^(^)-iE(^)(:^)^-^(^-^--)- (26) 

The summation over phonon modes yields 

A5 



J(A) 



2ATTn''pcl 

Substituting Eq.(27) into Eq.(19) we finally obtain 



(27) 



127rp;^ 



(28) 



Computing the spin matrix elements, one gets 



for the decoherence rate. Here we have introduced Ic = 
Ct/tOc- The generalization to finite temperatures is trivial. 
It results in the multiplication of the zero-temperature 
r of Eq.(28) by the factor coth(A/2r). Note that the 
form of Eq.(28) is independent of the form of the spin 
Hamiltonian as long as the crystal field is the only source 
of the spin-lattice coupling. The dependence of F on the 
parameters of the crystal field is absorbed into A = tkoc- 
The physics behind Eq.(28) is this The real-time 
quantum oscillations of S produce oscillating torque in 
the surrounding matter. That torque results in the tor- 
sion oscillations of the solid adjacent to the magnetic 
molecule, so that J = S + L = 0. The corresponding 
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time-dependent elastic deformation is a coherent quan- 
tum superposition of diverging and converging sound 
waves. That deformation is confined within distances 
of order l^. = Cf/cOc from S. By order of magnitude, Ic 
determines how far the sound goes away from the spin 
during the period of coherent quantum oscillations. At 
times of order the coherent superposition of sound 
waves is destroyed by a diverging wave of energy A that 
escapes to the surface of the solid. The energy UT de- 
termines the width of the excited state |1 >. It is easy 
to sec that TiT ^ L'^ /2Tc, where L = TiS is the angular 
momentum associated with the spin and 2c ~ plc is the 
moment of inertia of the volume of spatial dimensions Ic- 
If the magnetic molecule is to serve as a qubit, the 
frequency lVc should be sufficiently high, in the ballpark 
of 10^° — 10^^ s~^. Our study shows that it is possible 
only if the molecule is rigidly coupled to a solid matrix. 
This is contrary to the widely accepted view that the 
decoupling from the environment should be beneficial for 
the work of the qubit. In fact, any effort to make a loose 
connection of the spin with its atomic environment, by 
e.g. having a magnetic molecule on a some kind of a 
molecular leg, will most certainly decrease the tunneling 
spitting A. 

The rigid coupling to the atomic environment is also 

necessary to provide a small decoherence rate. Accord- 
ing to Eq.(28) the softening of the solid matrix with re- 
spect to shear deformations, that results, e.g., in the de- 
crease of Ci by a factor of four, increases T by a fac- 
tor of one thousand. For = 10 and Ct ~ lO^cm/s, one 
obtains from Eq.(28) r^lO-^s-i at Wc~10i°s-i and 
r^lO'^s"^ at Wc ~ 10"'^^s~^. These numbers show that 
magnetic entities can, in fact, be promising candidates 
for qubits. It should be emphasized that Eq.(28) estab- 
lishes the lower bound on the decoherence rate. Other 
effects unaccounted here, like interaction with nuclear 
spins, with free electrons, etc., can bring the decoherence 
rate up. 
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